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We present a microscopic study of the interlayer spacing d versus in-plane magnetic field B\\ 
phase diagram for bilayer quantum Hall (Qlf) pseudo-ferromagnets. In addition to the interlayer 
charge balanced commensurate and incommensurate states analyzed previously, we address the 
corresponding interlayer charge unbalanced "canted" QH states. We predict a large anomaly in the 
bilayer capacitance at the canting transition and the formation of dipole stripe domains with periods 
exceeding 1 micron in the canted state. 
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There now exists considerable experimental and 
theoretical [pj-p^ evidence for interlayer phase coherent 
states in bilayer quantum Hall (QH) systems at total 
Landau level filling fraction — ^- The ground state in 
these systems can be regarded as an easy-plane ferromag- 
nct 1^ , as a condensate of electrons in one layer and holes 
in the Landau level of the other layer jl^ , or as a super- 
fluid of Chern-Simons composite bosons [||. The compe- 
tition between tunneling energy, which pins the interlayer 
phase, and Coulomb interaction energy, which favors in- 
terlayer phase rigidity, yields a rich phenomenology. This 
is especially true when an in-plane magnetic field compo- 
nent B\\ , that favors the development of Aharonov-Bohm 
(AB) phases, is present. The nature of the ground state 
is dependent on the separation d between layers, on the 
strength of interlayer tunneling, which is conventionally 
parameterized by the splitting As^is between symmetric 
and antisymmetric bilayer single-particle states, and for 
AsAS ^ on B\\ . 

In this Letter we present a microscopic theory of the 
d — -Bji phase diagram, allowing for the spontaneously in- 
terlayer charge unbalanced "canted" commensurate (C) 
and incommensurate (I) QH states proposed recently by 
one of us ||l3|] , in addition to the corresponding charge 
balanced "planar" QH states. |^ Our main results are 
summarized by the phase diagram of Fig.^, where five dif- 
ferent phases occur. In the pseudospin ferromagnet lan- 
guage, the effect of tunneling is to add to the Hamiltonian 
an in-plane Zeeman pseudo-field, which winds uniformly 
in space along the x axis at rate Q = {B\\/ B±^{d/ for 
a field in the y direction. Here the magnetic length i 
is related to the perpendicular field by 27t£'^B± — hc/e. 
In the planar charge balanced commensurate (CP) QH 
state, the pseudospin magnetization faithfully follows the 
winding pseudo-field and its azimuthal angle is given by 
4>{r) — Qx. For sufficiently large Q, however, the cost 
in exchange energy of this variation becomes too large 
and phase-slip solitons are nucleated, leading to a incom- 



mensurate planar (IP) QH ground state, with the CP-IP 
transition in the universality class of the well-studied C-I 
transition . Within the Hartree-Fock (HF) theory we 
employ, the large d instability, previously studied only at 
Q = 0, is to a Wigner crystal state of the bilayer system 
1^,0. It has been argued previously Q that quantum 
melting of this crystal leads to a state which does not 
have broken translational symmetry but, like the Wigner 
crystal, is compressible, does not exhibit a QH effect, and 
is not distinguished by any symmetry from the state with 
uncorrelated v = 1/2 layers that is expected at very large 
d. We have followed that suggestion here by labeling 
this region "No-QHE" . Our work focuses on the small d 
instabilities of the planar QH pseudospin ferromagnets. 
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FIG. 1. Phase diagram illustrating the commensurate 
planar (CP), incommensurate planar (IP), commensurate 
canted (CC), and incommensurate canted (IC) QH states, 
for AsA5 = 0.126e^/e^. For smaller Asas, canted and in- 
commensurate states shift to smaller Q. The nature of the 
compressible "No-QHE" state cannot be established on the 
basis of HF theory. 
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At small d and intermediate in-plane fields, we find 
that drives the easy-plane XY anisotropy of these 
quantum Hall ferromagnets through zero, leading to a 
continuous Ising-like reentrant quantum transitions from 
the charge balanced CP and IP QH states to the corre- 
sponding commensurate (CC) and incommensurate (IC) 
interlayer charge unbalanced canted QH phases, recently 
predicted by one of us Our d — By phase diagram 

(Fig.|l|) was constructed by considering a set of single- 
Slater-determinant variational wavefunctions which allow 
translational symmetry to be broken in one direction: 

Mrhx]) = lliuxcy + vxcl^xM, (1) 

X 



where Cj, ^ and Cg ^ create electrons in top (T) 
and bottom (B) layers respectively, X is a Lan- 
dau gauge guiding center label in the lowest Lan- 
dau level, and the pseudospin magnetization rhx = 
(sin 9x cos (j)x , sin Ox sin (j)x , cos 9x ) at label X is the co- 
herent state label of the spinor {ux,vx), i-e. {ux, vx) = 
(cos[6'x/2],sin[6'x/2] exp[i(/)x]). In these wavefunctions 
4>x is the local pseudospin phase coherence angle and 6x 
is the local polar angle that specifies the magnitude of the 
charge imbalance between top and bottom layers. Tak- 
ing the expectation value of the microscopic Hamiltonian 
7i, this variational wavefunction leads to the microscopic 
energy functional Ehf[t^x] = 



EHF[r<Cix] = -^^SAS ^'""^X COs{<Px - Q^) + ^ [^^(^ " ' Ps{X - X')] COS Ox cosdx' 

X y x,x' 

Pd{X - X')smexsmex' cos{c^x - (l)X'), (2) 

y x,x' 



where l^sAS 



AO) 

^SAS 



yg^j-igg weakly with in- 



plane field, and we have assumed a presence of a fixed 
neutralizing background of positive charge. In the above, 
the terms proportional to H{X) arise from Hartree (elec- 
trostatic) contributions to the energy functional, while 
those proportional to Fs{X) and F]j{X) originate from 
Fock (exchange) interactions between electrons in same 
and different layers respectively. The first term in 
Ehf[i^x] is due to interlayer tunneling and incorporates 
the AB phases in the factor cos(0x — QX). For a model 
with arbitrarily narrow two-dimensional electron layers 

HiX)^ f^'-^^^^il^e^^-e-^'^'f^ (3a) 



2tt 



sm((l>x, 



2q 



Fc{X) 



~x^/2r 



dq 
2^ 



Vc{q,X/f)e' 



(3b) 



where VciqxiQy) ~ 2'iTe^/q and 2t: exp{~qd) / q for 
C — S and C = D respectively. Note that the exchange 
integral drops rapidly with orbit center separation, while 
the electrostatic integral falls only as X~'^ at large X, 
corresponding to interactions between lines of interlayer 
charge imbalance electric dipoles. 

The states we discuss are all extrema of this energy 
functional and therefore have pseudospin configurations 
that satisfy the following two equations: 



cos( 
cos 



H) 
9x) 



^SAS sin(QX) 
As AS cos{QX) 



T; Ex' Fd {X - X') sm{9x' ) sin(0x' ) 
t; Ex' Fd {X ~ X') sM9x' ) cos(0x' ) ' 
- X')-2H{X - X')] 



, cosi9x')[Fs{X 

M^x) ~ AsAS cos(,^x - QX) + -J- J2x' smiOx')FD{X - X') cos(0x 



(4a) 
(4b) 
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Eqs. ^ , ^ follow from the minimization of EHpliAix] 
with respect to (f>x and 9x, respectively. The numer- 
ator and denominator on the right hand side of Eq. ^ 
are the x and y components of the pseudospin Zeeman 
effective fields seen by the Hartree-Fock quasiparticles, 
which include contributions from interlayer exchange in- 
teractions. Note that the exchange local pseudo-field de- 
creases when (f>x is not constant. Equation 4b expresses 



the property that in the HF ground state the pseudospin 
is aligned at each X along the direction of the pseu- 
dospin Zeeman effective field. For example, the com- 



K,,iX,X') = — 



1 S^Ehf 



2iT Sm^^Sm^, 
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mensurate planar CP state {4>x = QX and sm{9x) = 1) 
solves these equations with a Zeeman pseudospin field 
AsAS + Ff){Q) lying in the xy-plane (with Fourier trans- 
form convention, f{p) = /(X) exp(— zpX)). Be- 
cause the same layer exchange energy normally domi- 
nates the electrostatic energy, the z component of the 
effective pseudo-field tends to have the same sign as the 
z component of the pseudospin, enabling the canted con- 
figurations we will discuss shortly. 

The CP and IP ground states are stable against cant- 
ing if all eigenvalues of the following matrix are positive: 

[2H{X-X')-FsiX~X')] 



2 



AsAs cos((/)x -QX) + ^J2^d{X- X') cos((^x - cj^x') 



X' 



(5) 



In the CP state, translational invariance simplifies the 
evaluation of the eigenvalue spectrum of K^ziX^X') 
which has plane waves exp(ipX) eigenfunctions with 
eigenvalues K^zip) = [Ag^s + Fd{Q) + 2H{p) - 
F5(p)]/(47r). An important feature of Kzz{p) is the non- 
analytic linear decrease in H{p) — {e'^d/2£'^){l — pd/2 + 
p^{(P — 3£^)/6 + . . .) at small p, which originates from 
the slow X~^ fall off in the (anti-ferroelectric) dipolc 
electrostatic interactions. One consequence is that the 
minimum in Kzz{p) always occurs at a finite p = p*, 
with p* (X (P at small d, (see the inset (a) of Fig.|^) lead- 
ing to a minimum at finite p. For a sufficiently large 
value of d, K* = Kzz{p*) is negative even at Q = 0; this 
critical value of d is assocated with the onset of the "No- 
QHE" compressible regime at Q = 0. For finite small 
Q, Fd{Q) = Fd{0) - AnpsiQi)^ + ... decreases, thereby 
expanding the "No QHE" regime quadratically with the 
applied in-planc field , as illustrated in Fig jl| 

The physics of the small d part of the phase diagram 
is different. The minimum of Kzz{p) {K*) occurs at a fi- 
nite, but much smaller value of p, and decreases with Q, 
crossing zero before the incommensurate state boundary 
is reached as illustrated in Fig. ^ The instability is asso- 
ciated with a change in sign of the pseudospin anisotropy 
energy and is closely analogous to a transition in which 
the easy axis of a thin-film ferromagnet changes from in- 
plane to perpendicular-to-plane, forcing the formation of 
stripe domains [l4| . In the quantum Hall case it is elec- 
tric rather than magnetic dipole interactions that force 
the transition to occur at finite wavevector. The CP- 
CC phase boundary is defined by K*{d,Q) = curve 
{d{Q) « (Sirpsi^ /e^){Q^ — ^~^), rising linearly above the 
d = critical value Qcp~cc — ?^^) as illustrated in 
Fig. ^ As indicated there, the CP-CC canting insta- 
bility is preempted at intermediate values of d and Q 
by the CP-IP commensurate-incommensurate transition, 
which for small Ag^s takes place at a small value of 
Q = Q CP-IP = 4/(77^), determined by the condition 
of vanishing soliton energy, with ~ y/inPps/ Asas 
the width of an isolated phase slip soliton in the in- 
commensurate state fisd - Within the HF approximation 
Qcp-cc IQcp-iP approaches 7r/4 as d ^ 0. At finite 
values of Ag^s, the CP-IP phase transition is located by 
extrema of the Hartree-Fock energy functional for which 
(j)x = (px— QX is periodic (modulo 27r) , at a Q for which 
the period extrapolates to infinity. 

Because of its inhoniogeneous nature, locating the 
canting instability of the IP state is significantly more in- 
volved, but conceptually similar to the analysis of the ho- 
mogeneous CP state discussed above. We first solve the 
mean-field equations for the IP extremum of the energy 



functional, explicitly seeking a self-consistent solution for 
(j)x with period a. This minimization is performed nu- 
merically for a finite value of Ly so that the number of 
distinct guiding centers per soliton Ng = aLy/ {2tt£'^) is 
finite. 

The stability limit of this planar soliton lattice QH 
state is marked by the appearance of a zero eigenvalue 
in a corresponding KzziX, X') matrix, which must be 
computed numerically. In analogy with the CP-CC tran- 
sition, here the transition is between planar (charge- 
balanced) IP and canted (charge-imbalanced) IC phases, 
both of which are incommensurate QH states, with IC 
state located on the small d side of this phase bound- 
ary, as illustrated in Fig. |l|. Some features of these non- 
trivial numerical results can be understood on the basis 
of qualitative considerations. ||l^ For large Q, the AB 
phases are so unfavorable for interlayer correlations that 
the ground state approaches a simple state in which tun- 
neling between the layers is ignored completely and (jfx 
approaches a constant to minimize the interaction energy, 
i.e., (px = QX. In this limit the canting stabilities must 
be the same as those of bilayer systems with As as — 0. 
Consequently, as illustrated in Fig. |l| we find that the 
canting transition is reentrant Jl^ , and for large in-plane 
fields only a single instability to a "No-QHE" state ex- 
ists, with the critical value of d universally smaller than 
that for the Q = Q case. The evolution of this IP-IC 
phase boundary with parallel field is calculated here for 
the first time, and it would be interesting to test exper- 
imentally. In the remainder of the paper we concentrate 
on the properties of the novel canted QH states, CC and 
IC, at small d, which we expect to be reliably rendered 
by the Hartree-Fock microscopic theory. 

Canted QH phases are distinguished from their pla- 
nar counterparts by a finite z-component of the pseudo- 
spin magnetization order parameter, ruz, which measures 
the interlayer charge imbalance nx — riB — TOz/(27r£^), 
that spontaneously develops inside canted states. As dis- 
cussed above, because the canting instability is at a finite 
wavevector p* (d) , for finite d the order parameter is 
staggered with period 2'k/p*. However, in the limit of 
small d i, such that p* — > (or looking at scales 
smaller than 27r/p*), is nearly uniform and, as in 
the magnetic case, it is often a good approximation to 
ignore the stripe domain structure and look at the p — 
(uniform 9x) extrema of the energy functional. Eq.pb] 
then leads to nonzero mz{d,B\\) — \/l — sin^ 9, with 
sine = AsAs/[Fs{0) - 2H{0) - Fd{Q)], which there- 
fore predicts the expected mean-field square-root growth 
of the canting order parameter (Fig.|^, inset (b)) inside 
the CC phase upon crossing the CP-CC phase bound- 
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ary in any direction. Inside the CC state, Eq.^ pre- 
dicts the quasi-particle gap (measured through the ac- 
tivated behavior of the longitudinal resistivity) to be 
Aqjj = As AS + Fd (Q) \/1 ^ (c^i ^11 ) 7 reduced relative 
to that of the CP state. 
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FIG. 2. Bilayer capacitance across CP-CC phase bound- 
ary, with peak's height set by wavevector p* (inset (a)) of 
pseudo- magnetization mz{X) that spontaneously develops at 
the transition. Inset (b) shows the increasing lowest harmonic 
amplitude of mz{X) and the corresponding wavevector p*{Q) 
that decreases inside the CC phase. 

We are now in the position to calculate the phase bound- 
ary for the CC-IC transition. Since the effect of cant- 
ing on the soliton physics is to reduce the one-particle 
tunneling energy by a factor of sin0 and the interlayer 
exchange interaction by a factor of sin^ 9, standard anal- 
ysis ll^] leads to the phase boundary Qcc-icid, B^^) = 

A/ (tt^o) /Vs'md, that is shifted to higher critical values 
of the in-plane field (see Fig.|l]). |l3) We also expect that 
because of the long-range (dipole) soliton interaction in 
the IC state, the usual l/|ln((5 — Qci\ rise in the 
soliton density will be replaced by a significantly slower 
\Q — Qcc-icl^^'^ increase inside the IC state. 

In contrast to the CC state, the charge imbalance rrix 
is a periodic function inside the IC state, even in the 
d — > limit, oscillating with period a of the soliton 
lattice around the mean value of the charge imbalance 
mg((i, 

The four QH phases that we have discussed are con- 
nected by novel continuous quantum phase transition dis- 
cussed in Ref. |l3j . We therefore expect and find a variety 
of universal experimental signatures near phase bound- 
aries of Fig.|l|. As illustrated in Fig.^ some of many strik- 
ing predictions is a strong d-dependent peak in the bi- 
layer capacitance near the CP-CC phase boundary, as 
well as the development of spontaneous interlayer charge 
imbalance (proportional to m^) inside the CC phase. We 



expect, that these and other critical anomalies as 
well as the dipolar stripe order of CC and IC phases, 
with period of order micron and tunable with d and | , 
should be readily observable in experiments that we hope 
our work will stimulate. 

In summary, we have computed d versus phase di- 
agram for QH bilayers, and found that in addition to the 
previously studied planar commensurate and incommen- 
surate phases (which at large d are unstable to a com- 
pressible "No-QHE" state), there exist interlayer charge 
imbalanced commensurate and incommensurate phases, 
in which pseudo-magnetization continuously cants out of 
the easy-xy-plane. We computed the charge imbalance, 
differential capacitance and single particle gap in these 
new phases, and suggested ways of accessing this physics 
experimentally. 
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